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Concordance between morphology and 
mitochondrial phylogenetic structure 
in Australian dobsonflies 
(Megaloptera: Corydalidae) 


A.M. Baker Et G. Theischinger' 

Abstract: A first-order comparison of morphology and mitochondrial phylogenetic structuring is conducted in re¬ 
presentative dobsonfly (Megaloptera: Corydalidae) taxa, occurring in the Sydney Water Supply Catchment, south¬ 
eastern Australia. The current morphological classifications of the taxa examined are concordant with relationships 
inferred from the phylogenetic data. However, estimated divergence times between pairs of taxa based on morpho¬ 
logy and biogeographic events, pre-date times inferred from genetic data. Possible reasons for this discrepancy are 
discussed, and the genetic results are compared with those obtained previously for other macroinvertebrate taxa oc¬ 
curring in the study area. 
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Introduction 

Freshwater catchments of south-eastern Australia 
possess generally rich and diverse macroinvertebrate fau¬ 
nas, although phylogenetic structuring within these as¬ 
semblages is, to a large extent, poorly known. 

The Sydney Water Supply Catchment (SWSC), in 
south-eastern New South Wales, encompasses the 
Hawkesbury-Nepean, Georges and Shoalhaven River 
Catchments and is governed by the Sydney Catchment 
Authority (SCA). While the primary function of the 
SCA is to supply water of sufficient quantity and quality 
to customers, it also aims to achieve the successful man¬ 
agement of genetic and species diversity, as well as main¬ 
tenance of prime habitat. In a previous study, BAKER et 
al. (in press 1) used genetic markers to describe the pat¬ 
tern and extent of phylogenetic structuring within repre¬ 
sentative macroinvertebrate genera in this region. Re¬ 
sults indicated a large amount of cryptic diversity within 
a range of crustaceans and insects, given the relatively 
small size of the study area. Distributions of most taxa ex¬ 
amined were associated with high altitude streams, above 
dam impoundments. These baseline data will be used in 
coming years to predict the likely consequences for bio¬ 


diversity in the SWSC, if certain developments go 
ahead, and provide information on how to moderate the 
impacts of such development. 

In the present study, we add to the previous work 
conducted on freshwater macroinvertebrates in the 
SWSC, by presenting preliminary analyses comparing 
morphology and mitochondrial genetic structure of cory- 
dalid megalopterans. To our knowledge, this study repre¬ 
sents the first research investigating phylogenetic rela¬ 
tionships among Australian megalopterans. 

Methods 

Taxa 

All Australian corydalids belong in the subfamily 
Chauliodinae. They are rather large insects, with both 
adults and larvae possessing strongly developed 
mandibles. The adults have long antennae, long legs, two 
pairs of functional membraneous wings and a very soft 
abdomen. The head of the larvae is flat and prognathous, 
the thorax is largely sclerotised, the legs are short, and 
there are lateral filaments (gills) on segments 1-8, with a 
pair of large prolegs on segment 10 of the abdomen. 


1 Unserem Kollegen und langjahrigen Freund (seit dem 1949 gleichzeitig von HA, in mehrfacher Damenbegleitung, und GT, solo, durchgeftihrten, his- 
torischen InsektennadehKauf bei Quinn Haslinger in Linz) Horst Aspock mit Respekt und in Anerkennung seiner wissenschaftlichen Leistungen ge- 
widmet. Andrew Baker und Gunther. Quirin Haslinger: Lehrmittel-Handlung mit einiger Tradition; man konnte dort von Sammlungsexemplaren und 
InsektenBiologien bis zu Insektennadeln und wissenschaftlicher Literatur alles Erdenkliche kaufen und bestellen, selbst Adressen einem unbekannter 
Kollegen bekommen. Eine schone Einrichtung, deren Auslagen Horst Aspock, ich (GT) und wohl auch andere „abnormale° Kinder mit viel Liebe und 
Interesse studiert haben. 
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Study Area 

The Sydney Water Supply Catchment is a network 
delivery* system drawing wafer from three major ouch' 
merits in south-east New South Wales, Australia: The 
Ha wkesbury/Nepean, Georges and Shoal haven. Howe- 
ver, the SWSC also encompasses several casrward-flo 
wing river catchments (Fig- 3), 


Sampling and Identification 

Representative larval samples of dobsonfly tax a from 
selected river catchments (Fig. 3), were obtained by 
kick'sampling riffle habitat. These samples were immer¬ 
sed under liquid nitrogen, returned to the laboratory and 
held in a - 80C freezer, until required for analysts. 




Fig. 1 : Archichoutiodes guttiferus (Walker), adult male 

Fig. 2: Protochautiodes biconicus Kimmiks, final instar larva, dorsal. 

Corydaltdae are generally associated with dear, cold 
water, usually in riffle situations, but at least some are 
known from swamps (Tj-iEISCHINGER 19*?I ). Adults (Fig. 
I) emerge in the warmer seasons, are short-lived, take lit* 
tie food and during the day are usually found resting on 
stems of plants bordering larval habitats. There is a 
strong crepuscular or nocturnal tendency of the adults, as 
well as an attraction to artificial lights. Copulation takes 
place on vegetation near the water, and egg'masses are 
laid on rocks or vegetation overhanging or projecting 
from water. After hatching, the larvae (Fig. 2) drop into 
the water where they live as predators under rocks and 
debris for 1-5 years. Pupation takes place on land in a 
simple chamber in soil or litter. The pupal stage lasts on¬ 
ly a few weeks. 



Group taxa were identified based on external mor¬ 
phology. using the available literature (TheisouNGER 
1999, 2000). with the aid of a high-powered stereo mi¬ 
croscope. 

Molecular Procedures 

Total genomic DNA was isolated using a modifica¬ 
tion of the CTAB/phenol-chloroform DNA extraction 
protocol (DOYLE & Doyle 1987). PCR was undertaken 
using mitochondrial coding region cytochrome c oxidase 
subunit I (COD primers ICO 1490 and HC02198 (Fol- 
MER et ah 1994). 

50 pt total volume PCR reactions contained: 55.8 
jjL ddH.O, 2 pL each of 10 pM primers, 1 pL 10 mM 
dNTPs, 2 jjL 50 mM MgCL, 5 pL I OX Reaction Buffer, 
0.2 pL DNA Tuq polymerase and 2 pL template DNA. 
PCR was performed on the Genoamp PCR System 9700 
(PE Applied Biosystems) with an initial 
denary ration step of 94 C for 3 minutes: 
30 cycles of 94 C for 30 seconds, SfrC tor 
30 seconds. 72’Cfor 30seconds; and a fi¬ 
nal extension step of 72 C for 5 minutes, 

PCR products were purified using the 
QIAquick PCR purification kit (Qiagen) 
and directly sequenced, using hah heavy 
and light strand primers, on an Applied 
Biosystems 377 automated sequencer, 
with complete overlap of sequences, to 
ensure DNA-stranJ homology. 

Sequences were aligned using BioEJit 
Version 5.0,9 (Hall 1999b Estimates of 
p-distance and net (corrected) pairwise 
sequence diversity ± SE (10,000 h*ot- 
straps) were calculated using MEGA Ver¬ 
sion 2.1 (Kumar ct a). 2001). Generic 
distances for calculation of net sequence 
diversity between taxa. were calculated 


Fig. 3: Sample sites [Key indicates site of 
origin for each haplotype). 
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using the DNA substitution model of best fit as determi- 
ned using Modeltest Version 3.06 (POSADA (Si CRANDALL 
1998). 

FELSENSTEIN’s (1988) likelihood ratio test, as imple- 
mented in the program TREE-PUZZLE Version 5.0 
(SCHMIDT et al. 2002), was used to test whether the COI 
data set adhered to a clock-like evolutionary model. 

As suggested by AviSE (1994), we used different phy¬ 
logenetic methods to ensure concordance among mit¬ 
ochondrial lineages; Maximum Likelihood (FELSENSTEIN 
(Si Churchill 1996), Maximum Parsimony (Eck & Da- 
YHOFF 1966) and Neighbor Joining (SAITOU (Si Nei 
1987). All phylogenies were inferred using PHYLIP Ver¬ 
sion 3.6a2 (FELSENSTEIN 1993). The Neighbor Joining 
method was invoked using the F84 model of DNA sub¬ 
stitution, gamma distributed rates, 5000 bootstrap pseu- 
doreplicates of the data (block size 3) and 10 jumbles. 
Block size was used to partition the sequence data into 
codons and the jumble option randomizes the input or¬ 
der of taxa during tree construction. The Maximum Par¬ 
simony method was invoked using Search for best tree, 
most thorough search, 1000 bootstrap pseudoreplicates of 
the data (block size 3) and 10 jumbles. The Maximum 
Likelihood method was invoked using empirical base fre¬ 
quencies, search for best tree, 1000 bootstrap pseudore¬ 
plicates of the data (block size 3), 10 jumbles of the da¬ 
ta, and three categories in the Hidden Markov Model. 

Results/Discussion 

Based on morphology, we found three species (two 
genera) of corydalids in our sample: Archichauliodes 
(Riekochauliodes) guttiferus (Walker) (A. guttiferus 
group), Archichauliodes ( Riekochauliodes ) plomleyi KlM- 
MINS (A. polypastus group) and Protochauliodes biconicus 
KlMMlNS. 

Identification of these three species is straightforward 
(Theischincer 1983, 1999, 2000). Adults of all 
Archichauliodes species (Fig. 4) can be distinguished from 
all Protochauliodes species (Fig. 8) by details of wing ve¬ 
nation; larvae, at least of the Australian forms, by partic¬ 
ulars of mandibular dentation and antennae (Figs 5, 9). 
The species groups of Archichauliodes ( Riekochauliodes ) 
differ from each other, particularly in the structure of 
adult terminalia and in the size of spiracular tubes on seg¬ 
ment 8 of the larvae (Fig. 6, 7). 



Fig. 4-6: Archichauliodes guttiferus (Walker): (4) adult, venation of 
forewing; (5, 6) larva, dorsal: (5) head; (6) abdominal segment 8. Fig. 7: 
Archichauliodes plomleyi group sp., larva, abdominal segment 8, dorsal. 
Fig. 8, 9: Protochauliodes biconicus Kimmins: (8) adult, venation of 
forewing; (9) larva, head, dorsal. 

Only representative COI gene sequences of each of 
the three taxa were analysed in this study. However, a 
comprehensive analysis of dobsonflies sampled from 
throughout the study area revealed no other corydalid 
taxa, with the majority of samples corresponding to 
Archichauliodes guttiferus. Also, the two species of 
Archichauliodes were sympatric at a number of sites in the 
study area (A. BAKER, unpublished data). 


Table 1. Polymorphic sites for each genetic type (Haplotype). Top row shows sequence position and a (.) indicates conformity to the reference sequence 
(Hi A. guttiferus). 

1 till! 1 1 » 1 1 I 1 1222 222222 22222 222223 33333 333333 3333 

1233445686891 1223456677889011 2344556866778B890001122334556677 

_ 381 03251 069161 4 4 6 8 4 02517368736 247956 14573 925870 36254 7395140328 

HI A guttiferus a g t a a a t t a t a t t t a a c c a t a t t a t t a t t t a a t a r a a t c t t c t c a t t t a a c a a t a t a a a t t t 

H2 A guttiferus .a. 

HI A guttiferus .a.c.gcc.c...c..t...taa.c.t..c.t.c..c....g...g...c.ctgt.cc 

H2 A plomleyi .a.t.g.c.c...c..t...taa.c.t.t.c..c....a...g...c.ct.t.cc 

H 1 A plomleyi taact at cacatgatta. aag.a . c.ctgccccg. ta.aatt .agt . a t t . tcc. a 
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Fig. 10: Midpoint-rooted Neighbor-Joining phylogeny. Scale indicates 
divergence and bootstrap percentages are printed at nodes. 

We describe five unique COI 486 base pair haplotypes 
(Table 1). There were no insertions or deletions detected 
in the sequences. Of the 486 base pairs of unambiguous 
sequence, 78 sites (16%) exhibited variation and 28 sites 
(6%) were phylogenetically informative. Of the 28 infor- 
mative sites, 1 (4%) was at the first position of the codon 
and 27 (96%) were at the third position. The transition 
(Ti) to transversion (Tv) substitution ratio was 0.8:1. 

Uncorrected pairwise p-distance sequence divergen¬ 
ces range from 0.3 - 1.1% within taxa and 3.8 - 14.2% 
between taxa (Table 2). Corrected net sequence diver¬ 
gences between pairs of taxa were as follows: A. guttife¬ 
rus - A. plomleyi (10.14% ± 2.86); A. guttiferus - P. bi¬ 
conicus (43.87% ± 13.11); A. plomleyi - P. biconicus 
(62.19% ± 16.10). 

In the present study, phytogenies were constructed 
under the F81 + G model of sequence evolution, with 
gamma shape parameter 0.0942. The Neighbor-Joining 
Phylogeny (Fig. 10) recovered three well-supported, di¬ 
vergent lineages of corydalid taxa. Maximum Likelihood 
and Maximum Parsimony phylogenies (not shown) reco¬ 
vered the same three lineages, with bootstrap support va¬ 
lues of at least 95%. Levels of mitochondrial divergence 
and the positioning of phylogenetic clades are concor¬ 
dant with expectations, based on current morphological 
taxonomic designations between Archichauliodes species 
and between Archichauliodes/Protochauliodes. 

The COI gene fragment did not exhibit any signifi¬ 
cant among-site rate heterogeneity (LRT (delta) = 4.26, 
df=3, P = 0.23). Gaunt & Miles (2002) calibrated the in¬ 
sect COI molecular clock, using the Blattaria/Orthoptera 
divergence, at 0.022% nucleotide substitutions per million 
years. Assuming this rate is representative of COI diver¬ 
gence rate in dobsonflies, and factoring the rate on the es- 


Table 2. Lower triangular matrix of uncorrected percent (p-distance) 
sequence divergence between pairs of haplotypes. 


timated net divergences, it can be predicted that: (1) the 
Archichauliodes species shared a common ancestor approxi¬ 
mately 5 million years ago, at the Miocene/Pliocene boun¬ 
dary, and (2) Archichauliodes guttiferus/plomleyi and Proto¬ 
chauliodes biconicus shared a common ancestor approxima¬ 
tely 24 million years ago, at the Oligocene/ Miocene 
boundary. Estimated levels of net divergence between A. 
guttiferus and A. plomleyi are the same order of magnitude 
as recorded levels of mitochondrial (COI) divergence 
among cryptic species of the freshwater shrimp Paratya (up 
to 8%), from the same geographic area (Baker et al. in 
press 2). Paratya and Archichauliodes are co-distributed 
throughout much of the SWSC, and the proposed Plioce¬ 
ne divergence within these taxa suggests that they are li¬ 
kely to have been affected by the same processes. 

Early Tertiary climate in Australia was warm and wet, 
but by the end of the Miocene Epoch the continent had 
become cooler and drier (CHRISTOPHEL & GREENWOOD, 
1989; POLE et al. 1993). By the Pliocene, the climate 
throughout much of inland Australia had become semi- 
arid, with continental aridity increasing during the Qua¬ 
ternary, punctuated by pluvial intervals (Alley 1998). Ba¬ 
sed on the molecular data, we may therefore propose an 
evolutionary scenario that ancestral populations of A. gut- 
tiferus/plomleyi and P. biconicus become fragmented associ¬ 
ated with ongoing aridity during the Miocene. Archichau¬ 
liodes populations may have subsequently become frag¬ 
mented during Pliocene aridity, evolved in isolation for a 
time, and recontacted during a Quaternary pluvial event. 

However, notwithstanding the concordance between 
currently accepted classifications and molecular results 
in the relative systematic position of the three studied ta¬ 
xa within Australia, the timing of divergence based on 
present-day geographic distributions of these genera, is 
inconsistent with the genetic data. According to the 
morphological classification and distribution as summari¬ 
zed by New & THEISCHINGER (1993), both Archichaulio¬ 
des and Protochauliodes occur in Australia and Chile, with 
Archichauliodes also in New Zealand, and Protochauliodes 
also in western North America. This suggests that a line¬ 
age including ancestral Archichauliodes was differentiated 
from a lineage including ancestral Protochauliodes , well 
before the break up of Gondwana, presumably still in the 
Mesozoic. This antiquity of the two genera based on their 
vicariant Gondwanic distribution, also suggests that the 
diversification of Archichauliodes into well-defined Aus¬ 
tralian species groups, occurred much earlier than indi¬ 
cated by the molecular results. 
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The use of molecular clocks in evolutionary biology 
remains controversial, for two major reasons. First, rates 
across lineages, genes and genomic regions may be 
heterogenous across taxa. Second, for taxa which diver- 
ged a long time ago, it is difficult to determine the actu¬ 
al number of nucleotide substitutions that have occurred 
in each lineage since common ancestry. This is because 
saturation (multiple substitutions at a single site) may 
mask the actual number of substitutions (ARBOGAST et 
al. 2002). Indeed, even under idealised conditions (i.e. 
with a clock specifically calibrated for the species under 
study using the fossil record), 95% confidence limits for 
divergence times estimated by a molecular clock, may 
still be quite large (HlLLIS et al. 1996). 

Two major factors may therefore have contributed to 
underestimation of divergence times, based on the pre¬ 
sent mitochondrial data: 

(1) the COl gene generally exhibits a relatively fast 
rate of evolution which may pronounce the effects of sa¬ 
turation, and 

(2) estimates were based on relatively few sequences 
from few taxa, which may reduce estimates of molecular 
diversity. 

Potential limitations of molecular clocks notwith¬ 
standing, a better understanding of evolutionary proces¬ 
ses affecting these genera could be attained by comparing 
representative sequences of all taxa within these two ge¬ 
nera, throughout their world-wide geographic range. 
This would indicate how closely allied the Australian co- 
rydalid taxa are with corydalids from South America, 
North America and New Zealand. In addition, sequen¬ 
cing of more slowly evolving nuclear and/or mitochon¬ 
drial genes in these taxa could reduce the effects of satu¬ 
ration and potentially allow the accuracy of genome-spe¬ 
cific calibration rates to be tested. 

In a previous study conducted in the SWSC, BAKER 
et al. (in press 1) reported high levels of genetic diversi¬ 
ty and cryptic species within atyid shrimp (Paratya) and 
leptophlebiid mayflies (Atalophlebia). Given the position 
of the corydalids as high-order stream predators, the rela¬ 
tively lower levels of diversity recovered in the present 
study are perhaps not surprising. 

However, it is not unusual for molecular, phylogenet¬ 
ic-based studies to inflate the number of recognised taxa 
in any given area. For example, over 100 new, old-world 
tree frog species from Sri Lanka were recognised after ex¬ 
tensive analyses of mitochondrial DNA, morphology, 
ecology and bioacoustics, where previously only 18 
species were known (MeeGaSKUMBURA et al. 2002). In 
this regard, and notwithstanding the observed concor¬ 
dance between morphology and mtDNA, it is important 
to recognise that the mitochondrial gene tree may be in- 
congruent with the organisma! species trees for the iden¬ 


tified dobsonfly lineages. Thus, the deep divergence and 
sympatry of mitochondrial lineages detected in the pres¬ 
ent study is certainly consistent with a scenario of his¬ 
torical genetic isolation, evolution along separate trajec¬ 
tories and subsequent recontact, but does not necessarily 
rule out post-separation (secondary) gene flow between 
them. Formal comparative analysis of nuclear (allozyme) 
markers would reveal whether diagnostic loci can be used 
as evidence for an absence of contemporary interbreed¬ 
ing, between sympatric lineages of corydalid taxa. 

The concordance between our molecular and mor¬ 
phological data as well as discordance between molecu¬ 
lar-based divergence times and biogeography, provide a 
fascinating starting point on the path to ultimately un¬ 
derstanding the processes which have affected the evolu¬ 
tion of this little known family of macroinvertebrates, in 
south-east Australia, or on Gondwana. It will be interes¬ 
ting to conduct a more extensive comparative analysis of 
evolutionary 7 relationships including other groups, parti¬ 
cularly Protochauliodes from Chile and Archichauliodes 
from Chile and New Zealand, in the future. It would al¬ 
so be valuable to conduct a cladistic analysis of morpho¬ 
logical data for the Chauliodinae, or better still, for the 
Corydalidae of the world. 
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